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RANDOM WALKS IN A DIRICHLET ENVIRONMENT 


NATHANAEL ENRIQUEZ AND CHRISTOPHE SABOT 


Abstract: This paper states a law of large numbers for a random walk in a random iid environment 
on Z d , where the environment follows some Dirichlet distribution. Moreover, we give explicit bounds for the 
asymptotic velocity of the process and also an asymptotic expansion of this velocity at low disorder. 

1. Introduction 

After a first breakthrough of Kalikow [3], giving a transience criterion for non-reversible multi¬ 
dimensional Random Walks in Random Environment, Sznitman and Zerner proved, several years 
later, a law of large numbers in [5], followed by a central limit theorem proved by Sznitman in 0. 
A generalization to the case of mixing environments was proved afterwards by Comets and Zeitouni 
in (we refer to m for an overview of the subject). Despite these progresses, many important 
questions, concerning recurrence or explicit criteria for a ballistic behavior, remain largely open. 

Among random walks in random environment, random walks in an iid Dirichlet environment take 
a special place, since their annealed law coincides with the law of some transition reinforced random 
walk having an affine reinforcement (see [2]). These reinforced walks are defined as follows. At 
time 0, we attribute, in a translation invariant way, a weight to each oriented edge of Z d , and each 
time the walk crosses an edge, the weight of this edge is increased by one. Finally, the walk is a 
nearest neighbour walk, which chooses, at each time, an outgoing edge with a probability which is 
proportionnal to its weight. 

The question of transience and recurrence for such walks, was answered by Keane and Rolles, 
in j2], in the case where the walk evolves on a graph which is a product of the integer line with a 
finite graph. In the context of trees, a correspondance between reinforced random walks and random 
walks in random environment was used before, by Pemantle, in |5J. Our purpose is to give some first 
results in the case of Z d . 

In this paper, we state a law of large numbers for such random walks, under a simple and explicit 
condition on the weights. Moreover, we give explicit bounds for the asymptotic velocity of these walks 
and also an asymptotic expansion of this velocity at low disorder. Low disorder corresponds, in the 
random environment model, to the case where the law of the transition probabilities is concentrated 
around its mean value, and, in its reinforcement interpretation, to the case where the initial weights 
of the transitions are large, so that these weights are not significantly affected during the life of the 
walk (at least, if the walk is transient). 

Let us precise that these walks do not enter the class of walks considered in and in several 

other works, asking the law of the environment to satisfy a uniform ellipticity condition. 

This ellipticity hypothesis is usually used in two ways : 

_ in the definition of Kalikow’s auxiliary Markov chain which involves the expection of the Green 
function of the walk killed when exiting a given set. The uniform ellipticity is then a comfortable 
assumption for checking the integrability of this Green function. 

_ in the estimates of the drift of Kalikow’s auxiliary Markov chain, the ellipticity condition often 
plays a key role. We overcome this difficulty by using an integration by part formula. 

In section 2, we give the definition of random walks in Dirichlet environment, remind their con¬ 
nection with transition reinforced random walks and we present our main results. In section 3, we 
present an integration by part formula that will be the key analytic tool in the proof of our results. 
Indeed, in section 4, it is shown how one can take advantage of the special form of the law of the 
environment, in order to estimate, using the formula of Section 3, the drift of the killed Kalikow’s 
auxiliary walk. 
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In section 5, we study the integrability of the Green function of the walk which ensures the 
existence of the original (non killed) Kalikow’s auxiliary walk and finish the proof of our first result 
by applying the law of large numbers of Sznitman and Zerner j!Jj. In section 6, we follow the scheme of 
|21 to get bounds for the asymptotic velocity of the walk, and deduce an expansion of the asymptotic 
velocity at low disorder. 


2. Definitions and statement of the results 


We denote by T 2d := {(xi, x 2d ) e]0,l] 2d , s.t., Ynti x i = 1}> and b Y {e-i)i<i< 2 d the family 
of unitary vectors of defined as follows: ( ei)i<i< d is the canonical basis of R d , and for all 
j £ }d T 1,..., 2d} ej — Cj— d . 

For all a := (ai, ...,a 2d ) £]0,+cx:[ 2d , we denote by A“ the Dirichlet probability measure on T 2d 
with parameters (an, ...,a 2 d) he. the measure on T 2d : 


r(ai + ...+a 2d ) a 1 

r(a 1 )...r(a 2d ) :Cl 


...x 


C*2d 1 
2d 


dx\...dx 2d -i. 


For a unit vector e of we will sometimes write, for reading conveniences, a e for the weight CQ 
where i is such that ej = e. 

Let us now introduce random walks in an iid Dirichlet environment on 7L d . 

We define an environment as an element ui = (w(x)) xeZ d where at any vertex x , u>(x) := (ui(x, x + 

ei), ...,u j(x,x + e 2d )) belongs to T 2d . We set /x := (8> A“, so that p, is a probability measure on the 

xez d 

environments such that (uj(x)) xe7 d are independent random variables of law A“. 

We denote by P“ the law of the Markov chain in the environment u> starting at 0 defined by: 

Vx € Z d , \/k € N, Vi = l,...,2d, P UJ (Xk+i = x + ei\Xk = x) = lu(x, x + ei). 

The law of the random walk in random environment (or the so-called annealed measure) is the 
probability measure P^ = J P u dfj,(cu). 

In we show that random walks in iid environment have the law of some reinforced random 
walk. The following proposition states that the case of a Dirichlet environment corresponds to a 
quite natural law of reinforcement: 


Proposition 1. The measure P^ satisfies that P^-almost everywhere, 

oti T Xi (h, X n ) 


P tJ ‘{X n+ 1 = x + e l |<r(A' fc , k < n)) = 


£fc=i + N k (n, X n ) 


where N(n,x) = {Nfin, x))i<j< 2 d and Ni(n,x ) = £" =0 l{ Xl +1 -x 1 =e i , x t =x}- 


We refer the reader to 0 for the proof. 


2.1. Bounds of the asymptotic velocity. We can now state our first result: 

Theorem 1. Let a := (or, ..., a 2d ) G]0, +oo[ 2d , and /i = (8) A“ a probability measure on the 

xez d 

environment. Let us assume that there exists i £ {1,..., 2d} such that a ei > 1 + ot- ei ■ 

The process X n is transient under P M , and 


3v € R d \ {0}, such thatP ^{—— —+ n ^oo v) = 1. 

n 


Moreover, for all i £ {1, ...,d}, 


rT-a o_ ei 1 

(£fef=i a fc) — 1 


< v.ei < 


d-ei Q—ei H - 1 

(Eli i«*)-i 


Remark 1: The assumption on ai ensures that the set 

d 

_J_ _|_ O!.— Gi 1, Oi— &i + l] 

2 = 1 

does not contain 0. It is a key ingredient in the check of Kalikow’s transience condition. 
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Remark 2: When the afi s are large, v becomes close to the vector —^- ^(a ei — a_ ei )ej. This is 

Efc—l ak i=l 

not surprising if one thinks at the corresponding reinforced walk: the initial weights of the transitions 
are large enough so that they are not significantly affected during the life of the walk, and the law of 

Qj. 

the walk becomes close to the law of the Markov chain with probability transition —- in the 

1 a k 

direction e^. 

Remark 3: In dimension 1, the condition of theorem 1 is actually optimal. Indeed from [J|, we know 
that the asymptotic velocity is not null if and only if either > 1 or > 1 , 

which corresponds exactly to a ei > 1 + a_ ei or a_ ei > 1 + a ei . Moreover, the asymptotic velocity 
of the walk is equal to — ,' a ~ ei _] . This shows the optimality of the lower bound in Theorem 1. 

2.2. Expansion of the velocity in the limit of large parameters. We turn now to the second 
result of the paper, which gives the asymptotic velocity of the walk in the limit of large parameters 
Ofc. Let us remind that, in the limit of large parameters ak, the environment is concentrated around 
its mean value. 

Let us fix some notations. We consider some fixed transition probabilities 

(mi) e T 2 d, 

and a parameter 7 > 0 (aimed to tend to 00 ). We consider the weights 

ak = jrrik, 

so that the expectation of the transition probability, E^(oj(x, x + e^)), is independent of 7 and equal 

to TTli. 

The mean environment (to^) defines the transition probabilities of an homogeneous walk on Z d , 
which is ballistic with asymptotic velocity 

2d 

dm = ' mfce/c, 

fc= 1 

when the mean drift d m is not null. We denote by G m its Green function. 

The following result gives an estimate in 0{\) of the asymptotic velocity (in section 6 we give 
explicit bounds for this estimate). 

Theorem 2 . Assume d m 7 ^ 0. 

For 7 large enough, Theorem 1 applies, i.e. there exists v 7 ^ 0 such that lim „_ >00 ^ = v , P M a.s.. 
Moreover, when 7 is large, we have the following expansion for v: 

v = dm- —(G m (0,0) - 1) + 0(4). 

7 7 

Remark 1: Surprisingly, the second order of the expansion is colinear to the mean drift d m - We see 
that (G m ( 0 , 0 ) — 1 ) > 0 , which means that there is a slowdown effect, since the second order term is 
directed in the opposite direction to the mean drift. 

Remark 2: In [Jj], the second author gave an expansion of the asymptotic velocity in the case of a 
uniformly elliptic environment. In this work, several of the estimates relied strongly on the ellipticity 
condition, so that the proofs of [S] have here to be modified. Nevertheless, if we apply the formula 
of to this case, we get the same expansion (many simplifications occur due to the particular 
expression of the covariance matrix). It is not surprising that the speed-up effect obtained in some 
cases of 0 is not observed in the case of a Dirichlet environment. The example of jji), section 2, 
was based, indeed, on some correlation between the transition probabilities in orthogonal directions. 
Here, there is a kind of independence of the transition probabilities in each direction, in the following 
sense: under p = ® A“, the law of cu(z, z + ef) is independent of the law of ( 1 4(71+4 

Remark 3: The Green function G m (0,0) has the following explicit Fourier expression 

1 f 1 
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(we refer to Step 2 of the proof of Proposition 3). 

3. An integration by part formula 


In this section, we present an integration by part formula on that will appear to be the key 
analytic tool in the estimation of the drift of Kalikow’s auxiliary walk. 

Lemma 1. For all a (E]0, +oo[ 2d , and all differentiable function f on 

2d 


fd\ a = 


'T 2d 


+ ... + ad 
Oil 


1 


I x\.fd\ a + — I x 1 .((y'x k -^-) - ^-)d\ a . 
'T 2d or J Tm ^ dx k dxi 


Proof: We use the well known identity between the Dirichlet law A“ and the law of the vector 
i •••> ) w ^ lere the random variables Zi are independent variables following the gamma 


distribution r(ai, 1) of density *e * on ®+- 

This identity implies 

f fdX s = 

JTo 


<T 2d 


/(; 


Z\ 


Z2d 


r(a 1 )...r(a 2 d) J Rf J Zi ’E?=i * 

Integrating by part with respect to z± , we get 

1 


)e ^i=i Zi z^ 1 dz\...dz 2 d- 


Ir 2d ^ r(ai + l)...r(a2d) 


f (/ — jr~) e ^= xZi zf x ...z^ 1 dz 1 ...dz 2 d 
J R 2 d dz i 


where f(z 1 ,...,z 2 d) ■= /(y^r 


E?Si 


Now, we decompose this last integral into the /-part and the part. 

Using, in the reverse sense, the “Gamma” interpretation of the Dirichlet law / \( ai+1, “ 2, ''' ,a:2 ' 1 ), the 
/-part becomes 

P(ai + ... + a 2 d + 1) 


T(ai + l)...r(a2d) Jr 2d J 1 

Now, the part writes 


f.xf 1 ...X 2 d d 1 dxi...dx 2 d-i = 


a i 


a 2 d 


Qi 


X\.fdX a . 


!T 2d 


air(ai)...r(a 2 d) 


f (zi.-^-)e 1 ■■■ z 2 d' i 1 dz 1 ...dz 2 d 

Jr i dzi 


and 


df = , Zl _ z\ 

^■g Zl ( E 2:d iZi (E ^ iZ . )2 


)/i - 


ZlZ2 


Z1 


2d 


(E-=i zd 2 


f2 - - - 


Z\Z2d 


(E?=i Zi) 2 


fd 


z■ 

2^i= i ^ 




Zk 


_^^2d _ 

k -1 2-,i=l 


where /*(*!,.., 2 2d ) = J^( ^ 2 f )• 

^ Ei=i ^ E»=i ^ 

The “Gamma” interpretation of the Dirichlet law (used for the third time) allows to 


conclude. 


□ 


4. Kalikow’s auxiliary walk 

We remind here the generalization of Kalikow’s auxiliary walk (see 0) which was already presented 
in 0. 

Let U be a connected subset of Z d , and 5 e]0,1]. We denote by dU the boundary set of U, i.e. 
dU := {z £ Z d \ U,3x £ U,\z - x\ = 1}. 
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For all z £ U and z' £ U U dU, and for all environment to, we introduce the Green function of the 
random walk under the environment w killed at rate 5 and at the boundary of U: 

G^ 5 (z,z') = E“ (fVl* fc=2 ,J 

\fc =0 / 

where Tjj = inf{k, Xk £ T, d \ U}. 

In the sequel, we will drop the subscript 6 when <5=1, and we will write Gf r (z,z') instead of 
Gu,i{z, z'). 

We introduce now the generalized Kalikow’s transition probabilities (originally, Kalikow’s transi¬ 
tion probabilities were introduced in the case 5 = 1): 


uu,s,z 0 {z,z + ef) = 


E u\ G u,, 5 (^ 0 , z)w(z, z + ef)} 
^n[^u,s( z 0 ’ 2 )] 


In order to give bounds for these transition probabilities, we will be led to apply the integra¬ 
tion by part formula of the previous section to the functions Gfj 5 (x,y), viewed as functions of the 
environment u>. 

For this purpose, we need the following lemma which gives the expression of the derivatives of 
these functions: 

Lemma 2. For all connected subset U of 1r, for all xi, £ 2 , £4 £ U, X 3 £ U U dU, |s 3 — X 2 \ = 1, and 
for all 5 €] 0 , 1 [, 

dGfj Ax 1 , X 4 ) 

Remark 1: The partial derivative is understood in the following sense: the function Gfj S can be 
defined by S n (Hu) n , where f Ijj is the transition matrix, defined in the proof below, whose entries 
are subjected to some stochasticity condition. But, at least locally, when 6 < 1, it can be extended by 
the same formula to a function of the variables ( flu(x,y )), which are not subjected to this relation. 
In this sense, the partial derivative has a clear meaning. 

Remark 2: When £3 £ dU , the right-hand term vanishes since Gfj 5 (x 3 , 0 : 4 ) = 0. 

Proof: Let us define the transition matrix Hu(x, y) = u>(x, y) if x £ U, and Qjj(x, y) = 0 if x £ dU. 
We have __ 

G^ 5 (x u x 4 ) = E< 5r Wf* 1 ,* 4) 

n >0 

and 

d ^ u \x u x 4 ) _ sx (o \ kl to 

d(w(X2,X3)) ^ U ( Xl '®2) u (x 3 ,x 4 ) 


so that, taking the derivatives term by term in the sum defining Gjf 5 (x 1 , 2 : 4 ), we obtain the 
result. □ 

We turn now to the estimation of the transition probabilities: 

Proposition 2. For all connected subset U of Z d , for all Zq,z G U, for all 5 G]0,1 [ and all i = 
1,..., 2d, 

2d 1 

■r /\ \ -1 7 1 ^ . Oii 

• */ (2_ ak) > 1, then j---- < uu,6,zo (z,z + ei )< ---- 

fc= 1 (z^fc= l a k) 1 \ 2 ^k=l ak ) 1 

2d 1 

• if C^oik) <1, then 0 < u> Ut g >20 (z, z + efj < / - 

fe=i (L fc =i ~ 1 

Proof: For the clarity of notations we give the proof for i = 1. 

Lemma 2 yields 

dGf, Azq, z) 

4 = 8 G% s (zo, z)Gfj 5 (z + *, z). 
d(u)(z, z + e*)) 
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We now apply Lemma 1 with / = Gf rs (zo, 2 ), viewed as a function of the only variables Xi := 
uj(z, z + ef) for i = 1, 2 d, and we get 


E n[ G uA z ^z)} = ai + "' + a 2 d E / tl G if,s(zo,zMz,z + e 1 )\ 


Oil 


1 

- E tl 


2 d 


Oil 

We recall that 


to(z, z + ei).Gjj^(zo, z) j <5 y ' ui(z, z + ek)Gfjg{z + ek, z) — dGfjg^z + ei, z ) 


fe=l 


(1) 


2 d 


6 u>(z, z + ek)Gij S {z + ek, z) — Gf r g(z, z) — 1 


fc=i 


so that the second term in the right side of ( 1 ) writes 


— Ep, [oj(z, z + ei).Gy g(zo, z) (Gy s (z, z) — 1 — SG^ s (z + ei, 2 ))] 

Oil 


so that we get 


E i*[ G u,s( z o,z)} = ^ + '^ + a 2 d E M [G%g(z 0 , z)w(z , z + ei)] 


OL i 


H- E^ [co(z, z + ei).Gy g(zo, z) (Gfj g(z, z) — 1 — 5Gfj g(z + ei, 2 ))] 

OL \ 


and for the ratio cbu,s,z 0 (z, z + ei) = 


£u,s,z 0 (z, z + e 1 ) = 


Oil 




E u 


Ep[u{z, z + eQG^fcp, 2 )] 

^[^(* 0 , ^)] 


;( 2,2 + ei).G£ 5 ( 2 0 , 2 ;) (G% iS (z, z) - <5G“ 5(2 + e u 2 )) 




Oj z )] 


( 2 ) 


But, 


2d 


u;(z, 2 + efc)(G^ 5 (2,2) — 5 Gfjg{z + ek, z)) — 1 


fc=1 


and therefore, for all k = 1, ..., 2 d, 


0 < Gu S (z, z) - SG^g(z + efc, 2 ) < — -——r. 

’ ’ uj(z,z + e 1 ) 


These inequalities allow to bound the ratio in the second term of the right side of (2), between 0 
and 1 , and this finishes the proof. □ 


5. Proof of Theorem 1 

We gather now all the ingredients of the proof of Theorem 1. We want to apply Sznitman and 
Zerner’s law of large numbers (Hj- From a careful reading of the proof of this law of large numbers, we 
can see that the only conditions that need to be fullfilled, are the integrability of the Green function 
Gy( 20 , 20 ) for all bounded U, and Kalikow’s condition. 

The integrability of the Green function is proved in the following lemma: 

Lemma 3. If there exists i £ {1, ...,2 d}, such that ai > 1, then for all connected subset U ofh d and 
all zq £ U, E /J ,[G‘i r (zo, 20 )] is finite. 

Proof: For the clarity of notations, we suppose that or > 1. 

Define now by N the least integer such that 20 + Nei belongs to dU. 

We have the following lower bound for the probability P(oj, Zq, U) to reach dU from zq without 
returning to 20 : 

IV-1 

P(v, zq, U) > ]^[ uj(z 0 + ke 1 , Zo + (k + l)ei). 
k =0 
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The number of returns to Zo before hitting dU, being a geometric variable whose parameter is 

precisely P{oj,zq,U), its expectation Gfj(zo,zo) is equal to —7 --. 

P(uj,zo,U) 


(nf=- 0 M*, + ke\, Zo + (k + l)ei)^ 


which 

□ 


We are now led to examine the integrability of 
is equal to (f T 37 d\ a ) N which is finite since a.\ > 1 . 

We now have to check Kalikow’s condition. 

We notice first that, under the assumption of Theorem 1, Lemma 3 applies and Kalikow’s auxiliary 
walk is well defined. Then, the monotone convergence theorem allows to make S converge to 1 in the 
inequalities of Proposition 2. 

We then deduce that the drift of Kalikow’s walk belongs to 

d 


1 


jj^ | [ttf:, Q—e-i 1 7 ^Gi (^ — ei T f] 


which does not contain 0, under the assumption of Theorem 1. This proves Kalikow’s transience 
condition. 

In order to estimate the asymptotic velocity of the process, we apply directly Proposition 3.2 of 
[H| which makes the link between v and the drift of Kalikow’s walk. 

Remark: in Lemma 3, we only got a sufficient condition for the integrability of the Green function 
to hold. A better result about this question would not have ameliorated the statement of Theorem 
1 as far as our check of Kalikow’s condition requires a stronger assumption. 

6. Proof of Theorem 2 


Theorem 2 of section 2 is actually a consequence of a more precise result, where the “O” in the 
expansion is replaced by an explicit upper bound. 

Let us fix some notations: we set 

2 d 

7 = 


i=l 


and 


O-i f / \\ 

mi = m ei = — = E (w(eij). 
7 


When 7 is large, the environment {co(x, et)) tends to concentrate around its mean ( rrii ), what can be 
seen from the expression of the correlations 


Cov M (w(a;,a; + ei),uj(x,x + ej)) = 


L , 


7+1 

i = j, 


7+1 


The mean environment ( 777 ) defines the transition probabilities of an homogeneous walk on Z d , and 
we define 


km — 2 'y ^ yjvn e 


m-e 


i=1 


so that 


1 ~ k m = y^(ym ei - A/m-eJ 2 


measures the non-symmetry of the walk. When km < 1, this walk is ballistic with asymptotic velocity 

2 d 

dm — ^ ^ 


i = 1 


and we denote by G m (-, •) its Green function. Let us define 


max^ 


T] m — 


1 kn 
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Proposition 3. Assume we are in the condition of application of Theorem 1, and that 

2 d 


7 


" I ]m 7 1; 


then we have the following estimate 


— dm(l-—r(G m ( 0 , 0 ) — 1 )) < 16 (-\ 

7 - 1 \7/ 


2 2 
V m 


7 / 1 _ 

' / 1 7 'In 


( 2 ) 


Proof: Considering the domain [/ = Z d , a killing parameter 6 < 1 and zq = 0, we get from formula 


us{z, z + ei) = mi + 


mi 1 E ll [G%(0,z)uj(z,z+ei)(G , g(z,z)-5G%(z + ei,z))\ 


7- 1 7- 1 Ep[G s { 0 ,z)] 

In the sequel, we will sometimes forget the superscript u> in G$, when there will be no ambiguity. 
Let us introduce a new probability on the environments fi(duj) given by 

OiM tl . 


We see that 

E„[G s (0, z)w(z, z + ei)(Gs(z, z) - 5Gg(z + e it z))\ 


= Ep,[(Gs(z, z) - SGs{z + ei,z))u)(z, z + e^)]. 


E »[Gs( 0 ,z)\ 

We proceed as in |H], and apply Kalikow’s formula (cf. the generalized version in {Bj, Proposition 
3.1) to the measure fi. 

It means that we have 

E[i[Gg (z, z)u>(z, z + ei)\ = Gg (z, z)w z {z, z + ef), 
where uj z is the auxiliary transition probability given by 

E iAGt{z,y)uj{y,y + ej)} 


*{y,y + ej) = 


E h[Gs(z,y)} 


Similarly, 


E p.[Gg {z + ei, z)lo(z, z + ei)\ — Gg (z + ei,z)ui“ e * (z,z + ei), 


where uj z+ei is the auxiliary transition probability given by 

&*+«(.. V + P )- E hl G s( z + e ^yMy^y + e j)} 

w \ y->y i tj )— r rn *)( i m 

E h[Gs{z + ei,y)\ 

Step 1 : We want to estimate the transition probabilities u z and ui z+ei . 
Lemma 2 yields 


d 


d 


duj(y,y + e k ) dw(y,y + ej) 


)G \5 (■, z) — 6G [5 (•, y)(Gg (y + e k , z ) - Gg(y + ej,z)), 


moreover 


2d 


J2u(y,y + e k ){G‘£{y,z)-5G l £{y + e k ,z)) = l y=z . 


k=1 


Using the integration by part formula given in Lemma 1, we get 


m ej E li [Gs(0,z)G s (z,y)} = E li [G s (0,z)Gg(z,y)uj(y,y + ej)] 

+~ E n [G«s(0,y) ( G s (y,z ) - 5Gg(y + ej,z) - l y = z ) Gg(z, y)u(y, y + ej)] 
7 

+~ E n [Gs{0,z)Gg(z,y) ( Gg(y,y ) - 6Gg{y + ej,y) - 1 )u(y,y + e j )\ 

7 
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But we have 
and if y 7 ^ z 


0 < v{y,y + ej)(G s (y,y) - 8G s (y + ej,y)) < 1, (3) 


\Gs(0,y)ui(y,y + ej)(Gg(y,z) - 5Gg(y + ej,z))\ < (2d - l)Gg(0, z ). (4) 

Indeed, for all k = 1, ■ • ■ , 2d, we have 

G$(v + e k ,z) (y,z), 

where T y is the hitting time of y (equal to infinity if the random walk never hits y). 
Since 

1-S£ k v{y,y + ek)EZ +ek [STy] = G * (2/ ’ y) ’ 

we get 


u(y,y + e k )(G%(y, z) - 6G%(y + e k , z)) < 


Gg(y,z) 

Gg(y,y)' 


But, we also have 


2d 


w ( 2 b V + e k)(Gg ( y , z ) - 5Gg (y + e k , z)) = 0 , if y ± z, 


k—1 


so that we have 


\^(y,y + ej)(G^(y,z)-SG^(y + ej,z))\ < ( 2 d-l) 


G$(y,z) 

Gfi(y,yY 


which immediately implies the estimate (4). 
The inequalities (3) and (4) imply that 


2 d. 


\m ej E^[Gs(0, z)Gs(z, y)\ - E tl [Gg( 0, z)Gg(z, y)w(y, y + ej)] \ < —E^[Gg( 0, z)Gg(z, y)]. 

This gives the following estimate for u> z 

I m ej -uj z (y,y + ej )\ < 

The same procedure gives the same estimate for A z+ei . 

Hence, we see that 

Eji[Gg(z, z)uj(z , 2; + ej)] = G™ +Am (z, z)(rm + A m(z, z + e*)), 
where A m(z, z+e ,;) is a correction to the homogeneous transition probability (mi) uniformly bounded 

by 

I Ami < —. 

7 

The same reasoning holds for 

Eji[Gg (z + ei, z)u>(z, z + ej)] = G™+ Am (z + ej, z)(n%i + A m(z, z + ej)), 
even if the correction term Ain is not the same. 

Step 2: We compare now the Green function c™+ Am with G™. 

This is done in 0, but we reproduce the main lines of the proof, since we want to obtain explicit 
bounds. We first introduce the symmetrizing function 

d 

r(z) = n 


i=1 


m- ei 


The Green function G™ is transformed into 


GT = M^G s Skm M^, (5) 
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where M$ is the operator of multiplication by <j>, and G s Sk is the Green function of the symmetric 
random walk with transition probability 

y/m ei m- ei 


Se„- — 5_ e . — 


2EttlVV-e 




with killing rates 5k m where 


2 d 


km = 2 E y/m^rn 


fe=l 


We refer to Step 2 of the proof of Lemma 4.3 in fBj for precisions about this fact. 
Hence, we see that 


1 


which means that 


We also have 


IG 


< 


1 — Sk n 
1 


■ 5fcll °° " 1 ~Sk m 


^m+Am _ ri 171 


- G™ = -G^ {I-(I- SAPrnGy)- 1 ) , 
where A P m is the matrix (A Pm)x,x+ ei = A m(x,x + e^) (and null anywhere else). 

Thus, we get 

Qm+Am _ G m _ SMy 1 GI^AP^y 1 G s Skm (/ - SM*AP m My 1 G s Skm ) M*, 

but 


\\M<j,AP m M± 1 || 00 < (max</>(ei)) — 

r i 7 


so that we get 


and 


(G7 +Am - G7)(x,y) < 4>{y - 


G7 +Am (^2/)<^(y-^) 


I — Zj* 1 

/ -L A/ m 1 — —//• 


1 — A* 1 7 Q 

1 i — —r/. 


This implies that, for all i £ {1, ...,2d}, 


u> s (z,z + ei ) -mi + ——-(G7(0, 0 ) -G?(ei, 0 ) - 1 ) 
7-1 


< 


1 


2 d 


1 


1 | 2 Vm 2 d ' - ‘Id 

7 “ 1 \ 7 1 - ™rjm 7 1 - frjr, 


2 d 
+ 2 — rj„ 


1 


< - 2 


% 


7 ^ I _ 2 d 
1 1 ~ '/n 


( 6 ) 


(we used here 17 ™ > 1 ). 
The sum 


2 d / \ 

E K - - 3 t(GT( 0 , 0 ) - G7(e,0) - 1) Vet 

i=i \ ^ ' 


^ 2d 

tends to d m (l-(G m (0,0) — 1)) when 5 tends to 1. Indeed, the sum > miG m (ei, 0).ei cancels, 

7—1 z —' 

' i=l 

due to the fact that for each i £ {1, ...,d} m ei G m (e,;,0) and TO_ ei G m (—e i; 0) are both equal to the 
common value y/m ei m- ei Gf.(ei, 0) (cf formula (5)). 

The triangular inequality combined with the 2d inequalities ( 6 ) gives that, for all z, 


limsup \\d& s (z) - d m ( 1--(G m (0,0) - 1))|| < 16 - - ^ 

1 7-! \7/ 1 — “ 
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where d^ 6 (z) = Y^t= 1 ^(z, z + ek)eu is the local drift of the transition probability uig. Proposition 
3.2 of 0 allows to conclude. □ 
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